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Abstract. For the stress analysis in a plastic body fi, we prove that 
there exists a maximal positive number C, the load capacity ratio, such 
that the body will not collapse under any external traction field t bounded 
by CYo, where Yq is the elastic limit. The load capacity ratio depends 
only on the geometry of the body and is given by 
1 J r \w\ dA 

7T = sup r \ ( \ i a t/ = HT-pII- 
C weLD(n)D J n \e(w)\ dV 

Here, LD(Q)d is the space of isochoric vector fields w for which the 
corresponding stretchings e(w) are assumed to be integrable and 70 
is the trace mapping assigning the boundary value 7x5 (w) to any w G 
I. !)[<>', -.. 



1. Introduction 

Consider a homogeneous isotropic elastic-perfectly plastic body Q. We prove 
in this study that there exists a maximal positive number C, to which we 
refer as the load capacity ratio, such that the body will not collapse under 
any external traction field t bounded by CY$, where Yq is the elastic limit. 
Thus, while the limit analysis factor of the theory of plasticity (e.g., P3E]) 
pertains to a specific distribution of external loading, the load capacity ratio 
is independent of the distribution of external loading and implies that no 
collapse will occur for any field t on dQ as long as 

ess sup \t(y)\ < CY . (1.1) 

yedQ 

Collapse will occur for some t if the bound is C'Yq with any C > C . 

The load capacity ratio depends only on the geometry of the body and we 
prove below that 

1 L \w\ dA 

C w eLD(n) D Jn\ £ {w)\ dV 

Here, LD{Q)d is the space of isochoric integrable vector fields w for which 
the corresponding stretchings e(w) are assumed to be integrable and 7/5 is 
the trace mapping assigning the boundary value ^d{w) to any w £ LD(Q)d- 
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The notion of load capacity ratio is an application to plasticity of ideas 
presented in our previous work @J E] where we consider stress fields on 
bodies whose maxima are the least. The general setting may be described 
as follows. 

Let Q represent the region occupied by the body in space so the body 
is supported on a part T of its boundary and let t be the external surface 
traction acting on the part T t of its boundary. Body forces may be included 
in the analysis using the same methods as in j^j but for the sake of simplicity 
we omit them here. The mechanical properties of the body are not specified, 
and so, there is a class of stress fields that satisfy the equilibrium conditions 
with the external loading. (Clearly, distinct distributions of the mechani- 
cal properties within the body will result in general distinct equilibrating 
stress distributions.) Each equilibrating stress field in this class has its own 
maximal value, and we denote by <7 opt the least maximum. 

Specifically, the magnitude of the stress field at a point is evaluated using 
a norm on the space of matrices. It is noted that yield conditions in plasticity 
usually use semi-norms on the space of stress matrices rather than norms. 
By the maximum of a stress field we mean the essential supremum over the 
body of its magnitude. Thus, we ignore excessive values on regions of zero 
volume. The traction fields that we admit are essentially bounded also. The 
set 17 is assumed to be open, bounded and its boundary is assumed to be 
smooth. Furthermore, it is assumed that T t and Tq are disjoint open subsets 
of the boundary whose closures cover the boundary, and that their closures 
intersect on a smooth curve. 

Subject to these assumptions (see further details in Section |2J our first 
result is 

Theorem 1.1. (?) The Existence of stresses. Given an essentially 
bounded traction field t onT t , there is a collection E t of essentially bounded 
symmetric tensor fields, interpreted physically as stress fields, that represent 
t in the form 



t-wdA= / a ij £ ij {w)dV, for all a € St, w £ C°°(17,M 3 ), (1.3) 



(ii) The Existence of optimal stress fields. There is a stress field a £ E t 
such that 





n 



where, e(w) = ^(Vu> + Vw T ). 




inf < ess sup \o-(x)\ 



ess sup \d{x) \ . 



(1.4) 



(Hi) The expression for <r° p . The optimum satisfies 
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where the magnitude of e(w)(x) is evaluated using the norm dual to the one 
used for the values of stresses. 

Item (i) above is of theoretical interest. It is a representation theorem 
for the virtual work performed by the traction field using tensor fields that 
we naturally interpret as stresses. It should be noted that the existence of 
stress is not assumed here a-priori. The expression for the representation by 
stresses turns out to be the principle of virtual work (| 1 . 3 1) . Thus, the equi- 
librium conditions are derived mathematically on the basis of quite general 
assumptions. Item (i) also ensures us that the representing stress fields are 
also essentially bounded. Item (ii) states that the optimal value is actually 
attainable for some stress field and not just as a limit process. 

Next, we consider generalized stress concentration factors for the given 
body. For a given external loading, traditional stress concentration factors 
are used by engineers to specify the ratio between the maximal stress in the 
body and the maximum nominal stress obtained using simplified formulas 
where various geometric irregularities are not taken into account. Regarding 
these nominal stresses as boundary traction fields, we formulate the notion 
of a stress concentration factor for a stress field a in equilibrium with the 
traction t mathematically as the ratio between the maximal stress and the 
maximum traction. Specifically, we set 

esssup xer2 |g(x)| 
esssup 3/ert \t{y)\ 

In particular, the optimal stress concentration factor for the given traction 
t, is 

opt 

K t = inf {K ti „} = ^ — - . (1.7) 

<res t esssupj, ert \t{y)\ 

Finally, realizing that engineers may be uncertain as to the nature of the 
external loading, we let the external loading vary and define the generalized 
stress concentration factor, a purely geometric property of the body Q, as 

K = su V {K t }, (1.8) 

t 

where t varies over all essentially bounded traction fields. In other words, K 
is the worst possible optimal stress concentration factor. Using the result on 
optimal stresses, we prove straightforwardly the following 

Theorem 1.2. The generalized stress concentration factor satisfies 

/ r \w\ dA 

K = sup r \ I \ I Air = H 7 °ll' ( L9 ) 
wec°°(n,R3) J a \ £ y w )\ av 



where 70 is the trace mapping for vector fields satisfying the boundary con- 
V 1 



ditions on Tn 1 



1 Further details on 70 are described in Section^ 
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To prove the theorems we use standard tools of analysis and the theory 
of LD-spaces given by [H1E3E1- Results analogous to Theorems Hl.l|) and 
H1.2|) were presented in our earlier work cited above. In [4], a weaker form of 
equilibrium is assumed, and in all earlier work we did not consider boundary 
conditions for the displacements on Tq. 

Next we turn to the adaptation needed for the application to plasticity 
It is assumed that the yield function is a norm on the space of matrices 
applied to the deviatoric component of the stress matrix. Thus, it turns 
out that the same mathematical structure applies if we consider isochoric 
(incompressible) vector fields in the suprema of Equations l|1.5p and (|1,9|) . 
For example, the analog of Equation 1|1.5[) is 



opt 

a t = sup 



J r t-wdA 



(1.10) 



where LD(Q)d is the collection of isochoric vector fields having integrable 
strains. 

It turns out that optimal stresses are related to limit analysis of plasticity. 
In fact, the limit analysis factor A* (see Remark 15.21 and P3E1IB]) is simply 
given by 

Furthermore, the expression for the optimal stress of Equation (|1.1()|) is im- 
plied mathematically by the results of Christiansen and Temam & Strang 
(PQE1IB])- This implies that the optimal stress fields do not require a special 
designs of non-homogeneous material properties but occur for the frequently 
used models of elastic-plastic bodies. In particular, elastic-plastic material 
will attain the optimal stress field independently of the distribution of the 
external load. 

We take advantage of these observations and introduce here the notion of 
load capacity ratio — a purely geometric property of the body. As described 
above, the load capacity ratio may be conceived as a universal limit design 
factor, which is independent of the distribution of the external loading. It 
immediately follows from its definition that 

C=l. (1.12) 

Section |21 introduces the notation, assumptions and some background ma- 
terial. In particular, the space LD(f2) of vector fields of integrable stretch- 
ings (or linear strains) (see (H1E3II1) is described. Following some prelimi- 
nary material concerning the boundary conditions in Section |31 the proof of 
the theorems is given in Section^ with some additional details in Appendix 
lAl The adaptation to plasticity theory, including the introduction of the 
load capacity ratio, is presented in Section [3 

It is noted that for structures, i.e., the finite dimensional approximations of 
the present setting, the expression for C may be set as a linear programming 
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problem. Thus, the load capacity ratio and the generalized stress concentra- 
tion factor may be computed approximately using standard algorithms. 

I would like to thank anonymous reviewers for pointing out the relation 
between the notion of optimal stress and the expression I obtained for it, and 
the limit analysis factor and the expressions for it in the work of Christiansen 
and Temam & Strang ([HEIIE])- 

2. Notation and Preliminaries 

2.1. Basic variables. We consider a body under a given configuration in 
space. The space is modelled simply by M 3 and the image of the body 
under the given configuration is the subset C I 3 . It is assumed that 
Q is open and bounded and that it has a C -boundary dQ. Furthermore, 
there are two open subsets Tq C dQ, and T t C dQ such that Yq is the 
region where the body is supported and Y t is the region where the body 
is not supported so that a surface traction field t may be exerted on the 
body on T t . Thus, it is natural to assume that T and Y t are nonempty and 
disjoint, To U T t = dQ, and A = 8Tq = dY t is a differentiate 1-dimensional 
submanifold of dQ. (The regularity assumptions may be generalized without 
affecting the validity of the constructions below.) 

Basic objects in the construction are spaces of generalized velocity fields. 
A generic generalized velocity field (alternatively, virtual velocity or virtual 
displacement) will be denoted by w. In the sequel we consider a number 
of Banach spaces containing generalized velocities and a generic space of 
generalized velocities will be a denoted by W. Generalized forces will be 
elements of the dual space W*. Thus, a generalized force F is a bounded 
linear functional F: W — > R, such that F(w) is interpreted as the virtual 
power (virtual work) performed by the force for the generalized velocity w. 
We recall that the dual norm of a linear functional F is defined as 

||F||=sup!^l. (2.1) 

i«6W IMI 

2.2. Virtual stretchings (linear strains) and stresses. As an exam- 
ple for the preceding paragraph, consider the case where W is the space 
L 1 (Q, R 6 ) of L 1 -symmetric tensor fields on Q. A typical element e £ L 1 (Q, R 6 ^ 
is interpreted as a virtual stretching field or a linear strain field. We will use 
|e(x)| to denote the norm of the matrix e{x). Various such norms are de- 
scribed in [3j. Thus, 

Hell! = J \e{x)\ dV. (2.2) 

Q 

The dual space L 1 (i7,M 6 )* = L°°(i?,R 6 ) contains symmetric essentially 
bounded tensor fields a that act on the stretching fields by 

o(e) = j a(x)(e(x))dV. (2.3) 
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Here, we use the same notation for the functional a and the essentially 
bounded tensor field representing it and we regard the matrix a(x) as a 
linear form on the space of matrices so a(x)(e(x)) = cr(x)ije(x)ji. Naturally, 
an element a £ L°°(i?,]R 6 ) is interpreted as a stress field. The dual norm of 
a stress field is given as 

= IMloo = ess sup \cr(x)\ . (2.4) 

x£f2 

Here, \cr(x)\ is calculated using the norm on the space of matrices which is 
dual to the one used for the evaluation of \s(x)\ (see PJ for details). Thus, the 
choice of the space L 1 (i?,M 6 ) for stretchings is natural when one is looking 
for the maximum of the stress tensor. 

2.3. The space of boundary velocity fields and boundary tractions. 

As another example to be used later, consider the space L 1 (r t ,M 3 ) of in- 
tegrable vector fields on the "free" part of the boundary. Its dual space is 

L\T t ,R 3 )* = L°°(T t ,R 3 ), (2.5) 

so that a generalized force in this case will be represented by an essentially 
bounded vector field t on T t . Using the same notation for the functional and 
the vector field representing it, we have 

t(u) = f t(y)-u(y)dA (2.6) 

so t may be interpreted as a traction field on T t as expected. The dual norm 
of the traction field t is 

\\t\\ = \\t\\oo = ess sup \t{y)\ , (2.7) 
y&t 

again, the relevant maximum. 

2.4. The space LD(f2) and its elementary properties. A central role in 
the subsequent analysis is played by the space LD{Q) containing vector fields 
of integrable stretchings (see (ZH3E3)- We summarize below its definition 
and basic relevant properties (see ^H] for proofs and details). 

2.4.1. Definition. For an integrable vector field w € L 1 (i7,IR 3 ), let Vw de- 
note its distributional gradient and consider the corresponding stretching (a 
tensor distribution) 

e(w) = ^(Vu> + V/). (2.8) 

The vector field w has an integrable stretching if the distribution e(w) is an 
integrable symmetric tensor field, i.e., it belongs to L 1 (i?,lR 6 ). For the sake 
of simplifying the notation, we use e for both the stretching mapping here 
and its value in the example above. The space LD(f2) is defined by 

LD{Q) = {w: Q -» M 3 ; w € L 1 (i?,lR 3 ), e(w) G L 1 (17, M 6 )} . (2.9) 
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A natural norm is provided by 

IMI = \\M\ld = IMIi + IkMlli (2-io) 

and it induces on LD{Q) a Banach space structure. Clearly, the stretching 
mapping 

e: LD{Q) — ► L 1 (i7, M 6 ) (2.11) 

is linear and continuous. 

2.4.2. Approximations. The space of restrictions to fi of smooth mappings in 
C°°(i7, M 3 ), is dense in LD{f2), so any LD-vector field may be approximated 
by restrictions of smooth vector fields defined on the closure Q. 

2.4.3. Trace mapping. There is a unique continuous and linear trace mapping 

7: LD(Q) — ► l}[dQ : M 3 ) (2.12) 

satisfying the consistency condition 

j(u\o)=u\ga (2.13) 

for any continuous mapping u £ C°(i?, M 3 ). Furthermore, the trace mapping 
is surjective. Thus, although LD-mappings are defined on the open set 17, 
they have meaningful L 1 boundary values. 

2.4.4. Equivalent norm. Let T be an open subset of df2 and for w € LD(f2) 
let 

\\w\\ T = J | 7 H| dA+||e(ti;)||i, (2.14) 
r 

then, ||iy||r is a norm on LD{Q) which is equivalent to the original norm 
defined in Equation 1)2.100 . 

3. Constructions Associated with the Boundary Conditions 

3.1. The space L 1 (9^,M 3 ) . Let L^dtt, M 3 ) C L x (df2,WP) be the vector 
space of vector fields on dfi such that for each u G L 1 (9l7,lR 3 )o, u(y) = 
for almost all y GTq. It is noted that the restriction mapping 

p : L x (dfl,'M?) — <-X 1 (ro,M 3 ), Po (u) = u\r (3.1) 

is linear and continuous. Thus, since 

L 1 (a^,M 3 )o = p 1 {0}, (3.2) 

L 1 (ai7,IR 3 )o is a closed subspace of L 1 (df2,R 3 ). 
The restriction mapping 

Pt : L 1 (a^,M 3 ) — ►£ 1 (T t ,R 3 ), |£>t(«)=«|r t (3.3) 

is also linear and continuous. In addition, as = <9r 4 = A have zero area 
measure, 

J \u\ dA = J \p t (u)\ dA, u E L 1 (9f?, M 3 )o, (3.4) 
an r t 
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so pt is a norm-preserving injection. 
Consider the zero extension mapping 

e : L\T t ,R 3 ) — L 1 (9f2, M 3 ) , (3.5) 

defined by 

eoW( y ) = | U(y) (3.6) 

V A 7 [0 fory^Tt. V 7 

Clearly, o eo is the identity on the space L 1 (r t ,IR 3 ). Moreover, for any 
u € L 1 (5l7,R 3 )o, eo(p t (u))(y) = u(y) almost everywhere (except for y € A), 
so eo o p t is the identity on L 1 ((9l7, R 3 )o- We conclude, 

Lemma 3.1. T/ie mappings pt and eo induce an isometric isomorphism of 
the spaces L 1 (9i7,M 3 )o and L^r^M 3 ). T/ie awa/ mappings e^ and p* induce 
an isometric isomorphism of the spaces L^I^M 3 )* and L l (dft, M 3 )q. Every 
element to £ L l (dft, M 3 )q is represented uniquely by an essentially bounded 
t e L°°(ri,IR 3 ) m ffce /orm 



t (u) = t-udA. (3.7) 



3.2. The space LD(ft)o of velocity fields satisfying the boundary 
conditions. Recalling the definition of the equivalent norm on LD(ft) in 
Equation (|2.14p . we set T = To in that equation. Henceforth, we will use on 
LD(ft) the equivalent norm 

HI = IMIr„ = / ItHI dA+ HeHlli. (3.8) 
r 

Consider the vector subspace LD(ft)o defined by 

LD(ft) = 7" 1 {L 1 (a^,M 3 ) } C LD(ft). (3.9) 

Thus, LD(ft)o is the subspace containing vector fields on ft whose bound- 
ary values vanish on To almost everywhere. Since 7 is continuous and 
L 1 (dft,M. 3 )o is a closed subspace of L l {dft, M 3 ), LD(ft)o is a closed sub- 
space of LD(ft). Combining this with Lemma (|3.1|) we obtain immediately 

Lemma 3.2. The mapping 

70 = Pt o i\ld{Q) ■ LD(ft) — L 1 (r t ,M 3 ) (3.10) 

is a linear and continuous surjection. Dually, 

7o = {l\LD{a )o T PV- L°°(r t ,R 3 ) — LD(^ (3.11) 
is a continuous injection. 

Observing Equation l|3.8|) . for each u> € LD(ft)o, 

\\w\\ = ||eH||i. (3.12) 
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Lemma 3.3. The mapping 

eo = e\ LD{Q)o : LD(n) -> L\f2, M 6 ) (3.13) 
is an isometric injection. 

Proof. Equation (|3,12p implies immediately that = ||e(w)||i for all w £ 
LD{Q)q. Being a linear isometry, the zero element is the only element that 
is mapped to zero, so £q is injective. In addition to relying on the technical 
property l|2.4.4jl of LD{Q) to show that £o is injective, it should be mentioned 
that this follows from the fact that for any vector field w on J?, e(w) = 
only if w is a rigid vector field, i.e., if w is of the form w(x) = a + b x x, 
a,b G M 3 . Now, the only rigid vector field that vanishes on the open set Tq 
is the zero vector field. □ 

4. The Mathematical Constructions 

Let t G L°°(r t ,]R 3 ) be a traction field on the free part of the boundary. 
Then, 7q(£) is an element of LD{Q)q representing t. The basic properties of 
elements of LD{Q)q are as follows. 

Lemma 4.1. Each S & LD(H)q may be represented by some non-unique 
tensor field a € L°°(i7,IR 6 ) in the form 

S = e* (a) or S(w) = J a(x)(e (w)(x)) dV. (4.1) 

Q 

The dual norm of S satisfies 

\\S\\ = inf || cr || oo = inf < ess sup \o~{x)\ > , (4.2) 

where the infimum is taken is taken over all tensor fields a, satisfying S = 
£g(<7), i.e, tensor fields representing S. There is a a G L°°{Q, R 6 ) for which 
the infimum is attained. 

Proof. The assertion follows from the fact that eo is a linear and isomet- 
ric injection as in Lemma Q3.3p and using the Hahn-Banach theorem. See 
Appendix El for the details of the technical lemma used and its proof. □ 

Applying this lemma to 5 = 7q(*) one may draw the following conclusions. 

Conclusion 4.2. Forces on the body given by essentially bounded surface 
tractions are represented by tensor fields on the body. These tensor fields 
are naturally interpreted as stress fields. The condition that a stress tensor 
field a represents the surface traction t is 

7 *(t)=e* (a), (4.3) 

and explicitly, 

J t- l0 (w)dA = J a(£ (w))dV, (4.4) 
r t j? 
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for each vector field w € LD(f2) , i.e., a vector field of integrable stretching 
satisfying the boundary condition on Tq. This condition is just the principle 
of virtual work which is a weak form of the equation of equilibrium and 
the corresponding boundary conditions. Thus, we have derived both the 
existence of stresses and the equilibrium conditions analytically under mild 
assumptions. 

It is noted that the subscript 0, only indicating the restriction of the vari- 
ous operations to fields satisfying the boundary conditions, may be omitted 
above. Also, as the restrictions of smooth vector fields on 17 are dense in 
LD(S1), it is sufficient to verify that the condition holds for smooth fields on 
Q. For such fields, the integrand on the left may be replaced simply by t ■ w. 

Conclusion 4.3. There is an optimal stress field a representing t and 



l|7o(*)ll = Halloo = inf < ess sup \a(x)\ > , (4.5) 
CT I xen ) 

where the infimum is taken over all stress fields a satisfying 7o(t) = £${o~), 
i.e., all stress fields in equilibrium with t. Thus, the infimum on the right is 
the optimal maximal stress cr° pt . In addition, by the definition of the dual 
norm we have 

hSWI = sup MWj^ (4.6) 

weLD(Q) \\ w \\ 

l*(7oM)l (A7 , 
= sup , (4.7) 

w eLD(Q) ¥{w)h 
where in the last line we used Equation 1)3.12(1 . We conclude that 



opt 

a t = sup 



f Tt t -70HCL4 



weLD(n) J n \e(w)\dV 



(4.8) 



Recalling that the restrictions of smooth mappings on Q are dense in LD(f2) 
and that for such mappings the trace mapping is just the restriction, the 
optimal stress may be evaluated as 



opt 

h = SU P 



f T t ■ w &A 



w 



fn\ £ ( w )\ dv ' _ (49) 

where the supremum is taken over all smooth mappings in C°°(1?,IR 3 ) that 
vanish on Tq. 

It is noted that the value of <r° p depends on the norm used for strain 
matrices. 



We now turn to the simple proof of Theorem 1)1.2 
Proof. We had 

opt l*(7oM)| I, * 

oy = sup = ||7 

wdLD{Q) \\ £ \ w )\\l 
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SO 



opt 

K= sup ^-r- = sup \ " '" Y" }> = ho* II = holl (4.10) 
<eL°°(r t ,R 3 ) ll r ll tei^Crt.R 3 ) I Pll J 

where the last equality is the standard equality between the norm of a map- 
ping and the norm of its dual (e.g., pp., 191-192]). □ 

5. Load Capacity for Plastic Bodies 

The analysis we presented in the previous sections may be applied to 
the limit analysis of plastic bodies. While in the preceding analysis the 
magnitude of the stress at a point was represented by the norm of the stress 
matrix, for the analysis of plasticity, the relevant quantity is the value of 
the yield function. The yield function is usually taken as a seminorm on the 
space of matrices — a norm on the deviatoric component of the stress. The 
necessary adaptation is as follows. 

5.1. Notation and Preliminaries. We denote by ttp the projection on the 
subspace of spherical matrices P, i.e., 

7rp(m) = ^mul, (5.1) 

and by ttd, the projection on the subspace of deviatoric matrices D so 

7T£)(m) = rri£) = m — 7Tp(m). (5-2) 

Thus, the pair (ttd,^p) makes an isomorphism of the space of matrices 
with D © P. We will therefore make the identification M 6 = D © P and 
M 6 * = D* © P*. We will use the same notation |-| for both the norm on M 6 , 
whose elements are interpreted as strain values, and the dual norm on R 6 *, 
whose elements are interpreted as stress values (although the norms may be 
different in general). Thus, we assume that the yield function is of the form 

Y(m) = \ir D (m)\. (5.3) 

For example, if we take |-| to be the 2-norm on M 6 * we get the Von-Misses 
yield criterion. In practical terms this means that the material yields at a 
point x when Y(a(x)) = Yq for some limiting yield stress value Yq € M + . 

Thus, we will extend the foregoing discussion to the case where Y ', evi- 
dently a seminorm, replaces the norm on the space of stress matrices. For 
the space of stress fields we will therefore have the seminorm || • ||y defined 
by 

\\ a \\ y = ||tt£) ° f ||oo = ess sup Y(a(x)). (5.4) 
The expression defining the optimal stress becomes 

cr° pt = inf ||erl|y = inf {||7r D o cr||oo} . (5.5) 

T(t)=s*(a) 7 *(i)= £ *(a) 
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The condition for collapse is <r opt ^ Y and we use ^ to denote the collapse 
manifold, i.e., 

* = {*| ^° Pt = ^0} • (5.6) 

Remark 5.1. The expression 

opt _ - n £ ||cr||y 
£ *(<r)= 7 *(t) 

for the optimal stress may be refomultated as follows. Recalling that Yq 
denotes the yield stress, we write a = o"i/A, ||<7i||y = Yq and noting that 
ll^/IMIHI^ = 1> we are looking for 

<r° pt = mf |ki/A||y (5.7) 



,cri£dB 

where B is the ball in L°°(J?, D) of radius Yo • Thus, 

a° pt = mf ^, (5.8) 

^ = su Pi A I e £ >i) = 7*(At)} • (5.9) 



opt 

a. F 



Clearly, in the last equation dB may be replaced by B because if we 
consider a with ||cr||y < 1, then, o\ = <r/||<7||y is in dB and the corresponding 
A will be multiplied by ||er||y < 1. 

In the last equation the unit ball B contains the (elastic) states of the 
material within the yield surface and we are looking for the largest multipli- 
cation of the force so the resulting stress is within the yield surface. Thus, 
we are looking for 

\ = \* = sup A, 3a £ B, e*(a) = 7* (At). (5.10) 



opt 

a 



•t 

which is the limit analysis factor (e.g., Christiansen [2] and Teman <fe 
Strang (HJ). 

The expression defining the generalized stress concentration factor as- 
sumes the form 

opt II II II I 

"t r P V r I'D <7||oo 

K = sup Tr-ii — = sup inf — — = SU P m ^ — fpn • (5-H) 

t ||*1| 00 t 7*(t)=E*(<r) ||t||oo t 7*(t)=e*(o-) ||t||oo 

For the application to plasticity, we use the term load capacity for C = 
l/K, Hence, 

1 . iitr 



C = ; opt/|m| =mf^. (5.12) 

For every loading t we set 
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so using <t?? = [|7*(A£)|| = Acr° p for any A > 0, one has 



opt 



Yn 



oo 



t\g ^U' opt " opt 
0V 0V 



IIMooAo- 



(5.14) 



't "i 

It follows that for any t, belongs to the collapse manifold \& and the 
operation above is a projection onto the collapse manifold. Thus, 



C = inf ■ 



(7 



opt 



inf ||i*||oo/*o, 



(5.15) 



and indeed, CYq = inft^ e * ||i*||oo is the largest radius of a ball containing 
only surface forces for which collapse does not occur. 

5.2. Constructions associated with the extension to plasticity. The 

top row of the following commutative diagram describes the various kine- 
matic mappings we used in the case of a norm on the space of stresses as 
considered in the previous sections. The subspace L 1 (J7, D) of L 1 (i7, M 6 ) 
contains deviatoric (isochoric) strain fields and there is a natural projection 
tt° d : L l (Q, R 6 ) -> L 1 ^,!)) given by ir° D (x) = tt d ° X- The inclusion of a 
subspace in a vector space will be generally denoted as i. We will also use 
the notation 

LD(f2) D =e 1 {L 1 (I?,D)} (5.16) 

for the subspace of isochoric L-D-vector fields. Thus, using Ed and jd for 
the restrictions £o\ld{Q) d and 1o\ld(Q) d ^ respectively, we have the following 
commutative diagram. 



L\r t 



The dual diagram is 
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ID 



LD(n) 



LD(Q) D 



£ " . T 1 



L L (f2, 



Lr{Q, D). 



(5.17) 



L°°(T t , 



L°°(T t , 



70 



Id 



LD{Q)l 



LD(Q)* D 



L°°(^,1R 6 ) 



L°°(f2,D). 



(5.18) 



Since Ed is just a restriction of Eq, it is still a linear, norm-preserving 
injection. Thus, the assertion of Lemma (|4.1|) and the subsequent conclusions 
hold where LD(Q)d, Ed, and 70 replace LD(f])o, eq, and 70, respectively. 
The expression for the optimal stress for the plasticity analysis is therefore, 



a. 



opt 



inf 

7 *(t)= £ *(<7) 



J r i • to d^4 



weLD(n) D 



(5.19) 
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Finally, the load capacity ratio is given by 

1 cr° pt f r \w\dA 

- = K= sup 1 ^ j -= sup T ^ 7 — — = || 7D ||. (5.20) 
^ teL°°(r t ,R 3 ) Plloo weLD{Q) D Jn\ £ \ w )\ dl/ 

Remark 5.2. Our result (|5.1i-)|l for the optimal stress associated with S = 
7 £ (*) ^ 



o-/ = sup I, -, ... . (5-21) 
w eLD(fi) D \\ £ {w)\\y 

For limit analysis in plasticity it is shown by Christiansen E] and Teman 
<fe Strang jH] that the kinematic version of the limit load is equivalent to the 
statical version above, specifically, 



A* = inf J sup {<■»(«,)} . (5.22) 
s(w)=i i aeB J 

We will show that the two expressions are equivalent for the setting of stress 
optimization. 

Equation 1)5. 21j) may be rewritten as 



°s Pt = SU P n -, \ M (5.23) 

w£LD(n)n,S(w)=l \\ £ K W )\\Y 

(5.24) 



mf 5(«;)=l Ikwlly 

where we used sup a ,(l/a;) = 1/ inf x. We conclude that 

1 



opt - *rf IK*- ( 5 - 25 ) 

cr 5 F S(w)=l 

On the other hand, if we replace in the kinematic version of the limit load 
the requirement a £ B by the requirement ||o*||y ^ Yq which is the analog 
in our setting, we obtain 

A* = inf J sup {e*{a)(w)} \ (5.26) 

S(w)=l [\\a\\ Y ^Y J 

= inf J sup {a(e(w))}\. (5.27) 

S(w)=l \\a\\ Y ^Y 



Using ||e(u>)|| = su P||o-||^i \ a { £ ( w ))\i we have 
so indeed A* = Yj/a^ 1 



X* = Y inf ||e(w)||y, (5.28) 

S(w)=l 
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Appendix A. The Technical Lemma 

Lemma A.l. Let W and V be two Banach spaces and tp: W — ► V an 

isometric injection. 

(i) For each S G W* there is some (non-unique) a G V*, such that 

S = <p*(a). (A.l 

(ii) The dual norm of S satisfies 

\\S\\ = inf \\a\\, (A.2 

where the infimum is taken over all a representing S, i.e., those satisfying 

S = ip*(a). (A.3 
(Hi) There is a a G V* such that 

\\S\\=mi\\a\\ = \\a\\. (A.4 

a 

Proof. Given S G W*, we may use the fact that 

ip^ 1 : Images C V — ► W (A. 5 

is a well denned linear isometry to write 

I-SO/?- 1 ^))! < iisihi^- 1 ^)!! = iishimi. (A.6 

It follows that 

S o ip' 1 : Images — ► R (A. 7 

is a bounded linear functional on the subspace linage^ C V. We recall that 
the Hahn-Banach theorem asserts that if U C V is a vector subspace and r 
is a bounded linear functional on U, then, r may be extended to a bounded 
linear functional a on V such that 

£7 (it) = t(u), for all u G U, (A. 8 

and 

■ r \cr(v)\ \t(u)\ , , . . „ 

|| CT || = SU pL^^ = SU pLA_^i = || T ||. (A. 9 

vev IMI «eu \\ u \\ 
Applying the Hahn-Banach theorem to the situation at hand, we conclude 
that the functional S o (p^ 1 may be extended to a linear functional a on V 
such that 

cr(u) = S o (^ _1 (u), for all u G linage^?, (A. 10 

or equivalently, 

S(w) =a(ip(w)). (A.ll 

By the definition of the dual mapping we conclude that S = f*(cr). 
In general, for any a G V* 

*/ mi W(f(w))\ \\a\\\\(p(w)\\ 

tp* (a) = sup ^ ^ Jl = sup 1 V ,;\| < sup " | "' \n ■ A -12 
*w IM| weW \\w\\ weW \\<p(w)\\ 

so 

\\cp*(a)\\ ^ \\a\\. (A.13 
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On the other hand, for any a G V*, such that S = </?*(cr) 

\S(w)\ \So^( v )\ \ a (v)\ 
sup ' y 1 = sup J 1 1 = sup 1 v , (A. 14) 

ingW ll^ll iiglmagei^ ||^|| i>£lmage</3 H^H 

so by the Hahn-Banach theorem 

||5|| = \\<p*(a)\\ = \\d\\, (A.15) 

where a is the element of V* extending S o <^ _1 and having the same norm 
as in Equation I) A .9(1 . 
We conclude that 

||5|| =inf \\a\\, (A.16) 

where infimum is taken over all a satisfying S = (p*(a). The infimum is 
attained for a as above. □ 
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